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ABSTRACT 

We report the detection of two planetary mass companions to the solar-type star HD 155358. The 
two planets have orbital periods of 195.0 and 530.3 days, with eccentricities of 0.11 and 0.18. The 
minimum masses for these planets are 0.89 and 0.50 Mj respectively. The orbits are close enough to 
each other, and the planets are sufficiently massive, that the planets are gravitationally interacting 
with each other, with their eccentricities and arguments of periastron varying with periods of 2300- 
2700 years. While large uncertainties remain in the orbital eccentricities, our orbital integration 
calculations indicate that our derived orbits would be dynamically stable for at least lO'* years. With 
a metallicity [Fe/H] of -0.68, HD 155358 is tied with the KlHI giant planet host star HD 47536 for the 
lowest metallicity of any planet host star yet found. Thus, a star with only 21% of the heavy-element 
content of our Sun was still able to form a system of at least two Jovian-mass planets and have their 
orbits evolve to semi-major axes of 0.6-1.2 AU. 

Subject headings: planetary systems — stars: individual (HD 155358) — techniques: radial velocities 



1. INTRODUCTION 

There is now strongly compelling evidence that plan- 
etary companions to main sequence stars are found 
prefere n tially around metal-rich stars. iGonzalej (|1997l 
Il998allbl . fl99l was the first to do a detailed comparison 
of the chemical abundances of stars with and without 
planets, and to conclude that there was a significant dif- 
ference in these samples. This finding has subsequently 
been supported by a large number of investigations, sev- 
eral of which included carefully selected control samples 
(e.g. [Re id 2002; Hcitcr & Luck 2003; Santos ct al. 2005.; 
Il'ischer fc Valentil 120051 ; iGrether fc Lineweaveii I2007D . 
Much of the early work in this area concentrated on 
determining the direction of causality in the observed 
correlation. That is, does the presence of planets pol- 
lute the stellar convective region causing an apparent in- 
crease in photospheric metallicity, or does the physics of 
planet formation significantly favor a metal-rich environ- 
ment? This issue has been addressed by searching for 
patterns of elemental enhancement in planet-host stellar 
photospheres that might result from the metal-rich (or H- 
and He-poor) detritus of planet formation. So far, there 
is little compelling evidence that post-formation stellar 
photospheric self-enrichment is responsible for the ob- 
served correlatio n; a primordial explanation is most likely 
(jGonzaled |2006[ ). Thus, searching for planets around 
stars with relatively low metallicity is extremely impor- 
tant for our understanding of the physics of planet forma- 
tion. We need to find the lowest metallicity stars around 
which planets are ab l e to fo rm. The radial-velocity sur- 
vey of iSozzetti et al.l (|2006l ) was the first to specifically 
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target low metallicity stars in order to understand in de- 
tail the \ow-Z end of the dependence of planetary system 
and formation on metallicity. 

We present here the discovery of two planetary com- 
panions to the star HD 155358. This is one of the sample 
of low- metallic ity stars that were i ncluded in our HET 
planet survey (jCochran et al.|[2004[ ) in order to begin to 
explore the low-metallicity tail of the heavy element dis- 
tribution of stars with planets. Details of the observa- 
tions are given in Section [2l Orbital solutions for the 
planets are given in Section 13.11 and 13.21 Section 13.31 de- 
scribes our use of a genetic algorithm to achieve a com- 
plete exploration of possible orbital solutions. The dy- 
namical stability of our preferred orbital solution is ex- 
plored in Section 13.41 Detailed parameters of the host 
star are given in Section [H and the implications of this 
fascinating system for planet formation theories are given 
in Section [51 

2. OBSERVATIONS 

All observations of HD 155358 were made in queue 
sched uled mode using the High Resolution Spectrom- 
eter (iTiliil [l998l) of the 9.2 m Hobby-Eberly Telescope 
(jRamsev et al.|[l998l) . A 400/zm optical fiber which sub- 
tends 2'.'0 on the sky tracks the star across the focal plane 
of the telescope. Starlight then passes through a molec- 
ular iodine absorption cell stabilized at a temperature of 
70C. The dense narrow I2 spectrum which is superim- 
posed on the stellar spectrum enab les us to model the 
instrumental point spread function fVale nti et al.|[T995f) 
and provides the velocity metric fo r pre cise radial ve loc- 
itv measurement (jButler et al.ll996l : lEndl et al.ll2000f) . A 
250/xm spectrograph entrance slit then gives a spectral 
resolving power of i? = X/5\ ~ 60,000. The spectrum 
is recorded on a mosaic of two 2048 x 4100 pixel E2V 
CCDs, which samples the spectrum at « 4 pixels per 
resolution element. The cross-disperser is set to obtain 
the spectrum between 4076A and 7838A. A small gap 
between the two CCDs falls at 5936A, allowing most of 
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the I2 absorption spectrum to be recorded on the "blue" 
CCD chip. 

We obtained a total of 71 high signal/noise radial ve- 
locity observations of HD 155358 between June 2001 and 
March 2007. All spectra were processed using standard 
IRAF routines, and velocities were computed using our 
high-precision radial velocity code. Table [1] gives the rel- 
ative radial velocities for HD 155358. The observation 
times and velocities have been corrected to the solar sys- 
tem barycenter. The uncertainty a for each velocity in 
the table is an internal error computed from the vari- 
ance about the mean of the velocities from each of the 
480 small chunks into which the spectrum is divided for 
the velocity computation. Thus, it represents the relative 
uncertainty of one velocity measurement with respect to 
the others, based on the quality and observing conditions 
of the spectrum. This uncertainty does not include other 
intrinsic stellar sources of uncertainty, nor any unidenti- 
fied sources of systematic errors. 

3. ORBITAL SOLUTION 

3.1. The First Planet 

A periodogram of the HET velocities of HD 155358 
show a very strong peak at ~195days, with sub stantial 
additional power in the 300-500 day range. The iScargld 
(fl982 ) false-alarm probability of this 195-day signal is 
4.5 X lO^^*'. We first fit a Keplerian orbit to the data 
us ing the GaussFit ge neralized least squares software 
of iJeffervs et al.l ()1988f ). This planet, with a period of 
193.8 d, eccentricity of 0.16, and K velocity of 31.0 ms~^ 
has a minimum mass of Msinz — 0.79 Mj and a semi- 
major axis of 0.63 AU. The RMS dispersion of the data 
around the orbital solution is 10.2 ms~^, and the reduced 
chi-squared = 9.67. This goodness-of-fit parameter is 
significantly higher than would be expected based, as it 
is based purel y on the internal errors quoted in Table [TJ 
iWrightl (|2005h gives a median "jitter" of 4.4 ms ^ for 
main-sequence stars similar to HD 155358. Wright also 
mentions that subgiants typically have jitters of about 
5ms-\ as do blue stars (with B-V < 0.6). HD 155358 
appears to fall into both categories (cf. section [H with 
our derived age of 10 Gyr and Hipparcos B-V=0.545. 
While some of the additional scatter might be explained 
as intrinsic stellar velocity variability expected in this 
star, there is probably some additional source for much 
of the scatter. Thus, we have searched carefully for ad- 
ditional planets in the system. 

3.2. The Second Planet 

A Lomb-Scargle periodogram of the residuals around 
the single-planet orbital fit shows very significant power 
over a broad range of periods from 200 to 800 days, with 
the strongest peak around 530 days. The false-alarm 
probability of this power at 530 days is less than 10^^. 
Smaller peaks around 330 days show FAP slightly less 
than 10"'^. Thus, we investigated the possibility of a 
second long-period planet HD 155358 c. We used Gauss- 
Fit to compute simultaneous least-squares solutions of 
double-Keplerian orbits to the observed velocities. We 
investigated periods around both 330 days and 530 days. 
While we were able to find formal Keplerian solutions for 
both periods, the solutions with Pc ^ 530 days gave sig- 
nificantly better fits than did shorter period second plan- 
ets. When we couple this with the results of our genetic 
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TABLE 2 
HD 155358 Double Planet 
Kepleman Solution 





Value 
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algorithm orbital fits (section 13. 3p and dynamical calcu- 
lations (section r3.4|) . we have concluded that the 530 day 
period is correct for HD 155358 c. This solution is given 
in Table [2] and is shown in Figure [TJ The uncertainties 
reported in Table [5] were generated by GaussFit from 
a maximum likelihood estimation that is an approxima- 
tion to a Bayesian maximum a posteriori estimator with 
a flat prior (Jefferys 1990). The observed HET velocity 
residuals to the HD 155358 c fit, phased to the period of 
HD 155358 b are shown in the upper panel of Figure [H 
and vice- versa in the lower panel. 

In fitting a Keplerian orbit to observed data, the ec- 
centricity is often the least well-determined of the orbital 
elements, and is the one element that is most often able to 
absorb any additional signals or periodicities that may be 
in the data. For example, in comparing the HD 155358 b 
single-planet orbit fit with the elements of HD 155358 b 
in the two-planet fits, the most significant changes in 
the orbital elements of HD 155358 b are that the eccen- 
tricity dropped somewhat and the K value increased. 
The uncertainty in the eccentricity of HD 155358 c is 
rather large. This is probably due to some regions of 
sparse phase coverage, particularly in the phase intervals 
0.30-0.45 and 0.8-1.0. The rms scatter of the observa- 
tions around this two-planet fit is 6.0 ms~^. This rms 
is quite consistent with a stellar internal "jitter" of 4.5 
to 5.0 ms~^ and our typical internal observed velocity 
uncertainty of 2-4 ms~^. If we add a stellar "jitter" of 
5ms^^ in quadrature with the internal errors given in 
Table [U we then get a reduced chi-squared for our two- 
planet fit of 1.15. 

3.3. Genetic Algorithm Investigation 

In order to understand fully the nature of the orbit 
of the second planet in this system, we used a standard 
genetic algorithm to explore further the x^-landscape of 
a 2 planet solution and to find possible additional 
minima. The initial orbital parameters were distributed 
randomly over a certain range of start values and during 
each iteration these parameters were randomly mutated. 
The solutions "evolve" by using as the environmental 
selection mechanism. Solutions that improve or at 
least do not degrade it by a certain threshold are allowed 
to survive and to multiply. Solutions that are above this 



threshold (e.g. < Xbcst + 0.1%), and are not able 
to move back into this range, are terminated after a 
few generations (typically 6 generations). Repeating this 
process many times (each time with new randomly se- 
lected starting values) allowed us to probe a much larger 
parameter space than is usually possible with standard 
orbital fit programs. 

For HD 155358 we performed > 10,000 trial runs, al- 
lowing the period of HD 155358 c to vary between 250 
and 2100 days (the time span of observations). We con- 
strained the period of HD 155358 b to the range of 185 
to 200 days and its eccentricity to eb < 0.5. 

Three minima were found by the genetic algorithm: 
the first with periods for HD 155358 c around 326 days, 
the second around with Pc 530 days and a third one with 
Pc around 1540 days. The parameters of the solutions 
found for the first two minima coincide with the solutions 
found by GaussFit. However, the third minimum at 1540 
days was intriguing, as it represents another possible pe- 
riod for HD 155358 c. A detailed examination of the so- 
lutions in this minimum revealed that their eccentricities 
were all in the range of 0.835 to 0.86. With such a high 
eccentricity the orbits of the two planets would cross and 
the system would quickly disrupt itself. It appears that 
this 1540 day solution represents a good formal fit for a 
physically impossible orbital configuration. In order to 
further rule out the 1540 days as the true period of the 
outer planet we performed 1000 additional trials in the 
period range of 1500 to 1600 days for HD 155358 c and 
using an upper limit on Cc of 0.6. This time no solutions 
was found with comparable low x^ values. We thus con- 
clude that the 530 day period is the most likely period 
for the second planet. 

3.4. Dynamical Stability Calculations 

While the orbits given in Table [2] represent the least- 
squares Keplerian orbit solution to the observations, they 
may not necessarily represent physically realistic orbits 
for two planets in this system. The inclusion of the sec- 
ond planet HD 155358 c drops the eccentricity of the first 
planet HD 155358 b shghtly, but the new outer planet so- 
lution has a somewhat non-zero eccentricity itself. While 
the orbits for this solution do not cross each other, they 
do have sufficiently close approaches that the planets may 
well be interacting dynamically. 

To investigate the orbital stability of the two-planet 
Keplerian orbital solution for the HD 155358 system, we 
conducted dynamical simulations using the N-body inte- 
grator SWIFTS A fuU description of SWIFT's capabil- 
ities is given bv lLevison fc DuncanI (|1994). We adopted 
the stellar mass M,=O.87M0 as discussed below in Sec- 
tion [H and the parameters of the two planets were those 
given in Table [2l The planets were assumed to be copla- 
nar, and the planet masses were taken to be the mini- 
mum values (i.e. sini = 1). All simulations which used 
shorter period in the range of 300 days discussed above 
in Section [3731 for the outer HD 155358 c planet resulted 
in ejection of one or both planets within 10^ yr. Thus, 
we reject all solutions with Pc in the 300 day range. The 
longer period (Pc — 530.3 days) for HD 155358 c is fa- 
vored dynamically, and was used in all subsequent simu- 

1 SWIFT is publicly available at 

|http: //www. boulder. swri.edu/~hal/swift.html| 



4 



Cochran et al. 



2002 



2004 



2006 



2008 



A 



> 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



HD 155358 




_L 



_L 



52000 



52500 



53000 53500 
JD - 2400000 



54000 



54500 



Fig. 1. — The HET data velocities of HD 155358 with the best double-planet fit (solid line). The RMS scatter of the data around the fit 
is 6.0ms~^. 



lations. 

The three-body system (star and two planets) was inte- 
grated for a total of 10^ yr. Initial orbits for the planets 
were taken from Table [H The system remained stable 
for the duration of the simulation. As shown in Figure[3l 
the two planets exchanged eccentricities on a relatively 
short timescale, with a period of approximately 2700 yr. 
The eccentricity of planet b varied between 0.02 and 0.18, 
while the eccentricity of planet c varied between 0.08 and 
0.22. The argument of periastron uj for both planets cir- 
culated with a period of about 2300 yr. The two planets 
spend most of their time with \ujh — luc\ between 120° 
and 240°. 

Noting the significant interaction between the two 
planets in the minimum-mass case, we tested the effect 
of inclination of the entire system by adjusting the plan- 
etary masses and repeating the simulations. More mas- 
sive planets (sin i < 1) should interact more strongly 
with each other, and thus may be less stable. Systems 
with sini > 0.33 (i >20°) remained stable for the du- 
ration of the tests (10^ yr). These results indicate that 
the dynamical stability evidenced in Figure [3] does not 
require the special circumstance of a nearly edge-on in- 
clination and thus minimum-mass planets. Even though 
the two planets are definitely interacting dynamically, 
the timescale of this interaction is long compared with 
the span of our observations. Thus, our two-planet least- 
squares Keplerian solution given in Section 13.21 is a valid 
approximation. 

Throughout the dynamical simulations, we have as- 



sumed that the two planets are in coplanar orbits. It 
is, however, conceivable that the planetary orbits are in- 
clined with respect to each other. We have conducted 
further tests with mutually inclined planetary orbits to 
determine the maximum mutual inclination for which the 
system remained stable. For the purposes of these exper- 
iments, the initial inclination of the inner planet was set 
to 0° and that of the outer planet was assigned a grid of 
starting values ranging from 5° up to 50° . The planetary 
masses were again assumed to be the minimum values, 
and the systems were integrated for 10^ yr. Systems with 
mutual inclinations less than 45° remained stable for the 
duration of the simulations. For higher values of the in- 
clination between the two planets, the eccentricities of 
both planets experienced chaotic variations resulting in 
the ejection of the outer planet within 10^ yr (45°) and 
2 X 10^ yr (50°). 

4. STELLAR PROPERTIES 

We measured the photospheric iron abundance, [Fe/H] , 
of HD 155358 by analyzing the "template" spectrum used 
in the radial velocity analysis. We assumed that [Fe/H] 
was an effective proxy for the general photospheric metal- 
licity, [M/ Hj. Our ana l ysis m ethod is the same as de- 
scrilDcd by iBean et all (|2006l ) for solar-type stars. We 
briefly describe the technique here, and refer the reader 
to that paper for a complete description. 

We fit synthetic spectra to the profiles of 30 Fe I lines 
in the observed spectrum. We generated the synthetic 
spectra with an updated version of the plane-parallel. 
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Fig. 2. — The best fit double-planet Keplerian orbit solution to the Hobby-Eberly Telescope data for HD 155358. The upper panel shows 
the phased inner planet (HD 155358 b) orbit, along with the observed data which have had the outer planet (HD 155358 c) orbit removed. 
The lower panel shows the phased HD 155358 c orbital solution, along with the observed data which have had the the HD 155358 b orbit 
removed. Two cycles are shown, with the second cycle in grey for clarity. The rms scatter of the data around this double-planet fit is 
6.0ms-l. 



local thermodynamic eq uilibrium (LT E), stellar analysis 
computer code MOOG ljSnedenlll973D . We assumed as- 
trophysical log gf values for the Fe I lines, which were 
determined by fitting a solar spectrum and assuming the 
solar iron abundance log e(Fe)0 — 7.45. Our analysis 
was therefore differential to the sun. 

We also adopted model atmospheres computed with 
the general-p urpose stellar atmosphe re code PHOENIX 
(version 13, iHauschildt et al.l Il999f ) for our analysis. 
The model atmosphere effective temperature, Te//, was 
constrained using the {{B — V), [Fe/H ]| - T^jf rela- 
tionship of iRamirez fc Melended ()2005( ) . The surface 
gravity, log g, was constrained by interpolating the 
IBertelli et al.l ()1994f ) evolutionary isochrones for a given 
My, Teff, and [Fe/H]. We took the needed V magnitude, 

1 B — V) colo r, and parallax from the Hipparcos catalog 
Perrvman et al.lll99 7). 
We determined [Fe/H], microturbulence ^, and macro- 



turbulence 1] for HD 155358 by fitting the line profiles 
in the observed spectrum simultaneously. As the con- 
straints on the stellar Tej / and log g described above are 
both dependent on [Fe/H], these parameters were also 
varied simultaneously. We fit the observed spectrum by 
using an a daptation of the "Marquardt" mi nimization 
algorithm ()Marquardtlll963l : iPress et al.lll986[ ). The un- 
certainty in the determined [Fe/H] was calculated from 
the scatter in the abundances determined for each line 
individually and deviations due to the uncertainties in 
the adopted T^f / and log g. 

In addition to the standard spectroscopic parameters, 
we have also estimat ed the mass and age of HD 155358 
by interpolating the IBertelli et al.l (|1994f l evolutionary 
isochrones as we did to constrain the log g. The pa- 
rameters that we determined for HD 155358 are given in 
Table |3l Table |4] lists the literature values of the stel- 
lar parameters for comparison. Most notably, we found 
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Fig. 3. — Results of dynamical simulations of the HD 155358 system. The evolution of scmimajor axis (a), eccentricity, and the argument 
of periastron (uj) are shown for a typical 10'' year interval. The simulation was run for a total of 10** years. The two planets exchange 
eccentricities on a timescale of about 2700 yr, and lu circulates with a period of about 2300 yr. The semimajor axes of both planets are 
stable throughout the simulation. These computations were for a star of mass 0.87 M©, and planet masses fixed at the minimum mass 
(sini = 1). However, similar results arc obtained for sini > 0.4. 



[Fc/H] = -0.68 ± 0.07. This is in very good agreement 
with the literature values, which have a range -0.72 < 
[Fe/H] < -0.60, and median value [Fe/H] = -0.67. The 
derived value of \ogg, coupled with the mass, age, and 
effective temperature of HD 155358 indicate that the star 
has evolved off the main sequence. 

5. DISCUSSION 

The planetary system around HD 155358 is quite re- 
markable in that it contains two Jovian-mass planets 

in relatively short-period orbits around one of the most 
metal-poor planet-host star yet found. Thus, this system 
is crucial to our understanding of the physic;s of planetary 
system formation and early dynamical evolution. 

The final configuration of any planetary system results 
from a competition of several different physical processes, 
each occurring on its own independent timescale. The 
planets themselves must form from the remnant disk of 



TABLE 3 
Derived Stellar Parameters 
FOR ilD ir,r,:]r,s 



Parameter Value 



Teff 5760±101 K 

log g 4.09±0.10 (cgs) 

[Fe/H] -0.68±0.07 

(, 1.55±0.15 km s"' 

7) 1.89±0.20 km s"! 

M* 0.87±0.07 Mq 
Age ~ 10 Gyr 



material left over from star formation. This must happen 
before the circumstellar disk is dissipated. After planet 
formation, the system must evolve into the final state 
in which we detect it. The dependence of each of these 
processes on the metallicity of the material from which 
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TABLE 4 

Stellar, Parameters for HD 155358 from 
Other Sources 



Ti,ff log g [Fe/H] Mass Age Source 
(K) (cgs) (Mq) (Gyr) 



5760 


4.09 


-0.68 


0.87 


10.1 


1 


5870 


4.19 


-0.67 


0.75 




2 


5868 


4.19 


-0.67 






3 


5914 


4.09 


-0.61 






4 








0.84 


10.2 


5 


5860 


4.10 


-0.60 






6 


5818 


4.09 


-0.67 


0.89 


12.7 


7 


5888 




-0.69 


0.83 


13.8 


8 
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5808 




-0.72 


0.89 


11.9 
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the entire system forms then determines the range of 
al lowed planeta r y syst ems as a function of metallicity. 

iMatsuo et alJ ()2007t ) investigated the circumstances 
under which stars of various metallicity and mass could 
form planets within the framework of both the core- 
accretion model and the disk-instability model. They 
derive a low metallicity threshold for planet formation 
within the core accretion model of [Fe/H] = -0.85 for a 
disk of 5 times the mass of the minimum mass solar neb- 
ula (MMSN) and [Fe/H] = -1.17 for a disk of 10 times 
the MMSN. They also showed that the disk instabihty 
model has essentially no dependence on stellar metallic- 
ity, but would require a disk mass of at least 10 times the 
MMSN to form these planets. Thus, the two giant plan- 
ets we have discovered in the HD 155358 system could 
have formed from either formation mechanism. In either 
model they require a disk which was substantially more 
massive than the minimu m-mass solar nebula. 

Ilda k LinI ()2004l . l2005l ) investigated the dependence of 
the formation of planetary systems on stellar metallic- 
ity within the core-accretion model. The motivation for 
their model was to explain the observed correlation be- 
tween high stellar metallicity and planetary systems eas- 
ily detected by RV techniques. The basic idea is that 
in systems of high metallicity, the dust surface density 
in the protoplanetary disk will be enhanced. This will 
lead to pla netary core s being formed on a much shorter 
timescale. Ilda fc LinI f|2005) derived a planetary core 
mass at time t (prior to depletion of the feeding zone) 
that is proportional to lO'^'^, where Z is the stellar metal- 
licity (1O[^°/^1*Z0). For a giant planet to grow, a core 
must achieve a mass above a critical value of about 10- 
2OM0 before the gas disk is dissipated. If the core does 
reach a super-critical mass while there is still significant 
gas content in the disk, it will rapidly accrete gas and 
grow to become a gas-giant planet. If the core remains 
sub-critical through gas disk dissipation, then the plane- 
tary core will remain a Neptune or a super-Earth. Since 
the disk dissipation has a much weaker (if any) depen- 
dence on stellar metallicity, gas-giant planet formation 



around low Z stars is significantly more difficult than 
around stars of solar or higher metallicity. Unless the 
circumstellar disk is rather massive to begin with, there 
simply isn't time to grow critical-mass cores before the 
gas disk dissipates. Thus, one might expect a critical 
lower metallicity threshold for the formation of gas giant 
planets. The value of this threshold would depend on the 
physics governing the overall mass of the disk, which are 
poorly understood. 

The third relevant timescale for determining the final 
configuration of a planetary system is that of the plan- 
etary system dynamical evolution, which transforms the 
system from the one which formed the planets into the 
system which we detect today. There are two basic phys- 
ical processes governing post-formation dynamical evolu- 
tion. These are tidal migration prior to disk dissipation, 
and planet-planet or planet-planetessimal scattering fol- 
lowing disk dissipation. Only tidal migration can have 
any significant dependence on stellar metallicity. The 
disk evolution is driven by the disk viscosity. In type 
II migration, in which the planet has opened a gap in 
the disk as a result of depleting the local gas in the disk 
from runaway gas accretion, the planet is l ocked into 
the viscous evolution of the disk. Ilda fc LinI (|2005[ ) give 
a disk migration timescale that depends linearly on the 
disk viscosity, a. The dependence of a on metallicity 
is not known, and most models assume a constant a of 
order 10"''^ for all disks. Thus, it appears that once a 
system of several planets is able to form in a disk, the 
subsequent dynamical evolution of that system is prob- 
ably reasonably independent of Z. Detailed modeling of 
the dependence of disk evolution and type II migration 
on Z is needed. If the dynamical evolution of the plan- 
etary system is due to tidal migration in the disk, this 
migration must still occur before the disk is dissipated. 
If the final configuration of the HD 155358 system is due 
to post-formation planet-planet scattering, or dynamical 
interactions of the planets with remnant planetessimals, 
then we would not necessarily expect any significant de- 
pendence of the process on the stellar metallicity. For 
low metallicity systems, the theoretical challenge is to 
form the planets of the required mass before the disk 
dissipates. 



This material is based on work supported by the 
National Aeronautics and Space Administration under 
Grants NNG04G141G and NNG05G107G issued through 
the Terrestrial Planet Finder Foundation Science pro- 
gram. The dynamical simulations used in this work were 
carried out at the Texas Advanced Computing Center 
(TACC). The Hobby-Eberly Telescope (HET) is a joint 
project of the University of Texas at Austin, the Penn- 
sylvania State University, Stanford University, Ludwig- 
Maximilians-Universitat Miinchen, and Georg-August- 
Universitat Gottingen. The HET is named in honor of 
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